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Математичне моделювання процесів, 
які мають розриви, за допомогою розривних 
інтерполяційних сплайнів 


У роботі запропонований метод побудови розривного інтерполяційного лінійного сплайна для наближення 
функції однієї змінної, що має хоча б один розрив першого роду у вузлах розбиття області визначення 
функції, використовуючи метод мінімакса. Визначений загальний вигляд похибки наближення функції 
побудованою розривною конструкцією в інтегральному вигляді, та наведені оцінки похибки наближення в 
кожному інтервалі розбиття. Запропонований метод можна буде використати для математичного моделю- 
вання розривних процесів в медичних, геологічних, космічних та інших дослідженнях. 


Вступ 


Задачі дослідження процесів, що мають розриви, виникають значно частіше, 
ніж задачі дослідження неперервних процесів |1|, 12). Наприклад, при дослідженні 
внутрішньої структури тіла корисно враховувати його неоднорідність, тобто різну 
щільність в різних частинах тіла (кістки, серце, шлунок, печінка тощо мають різну 
щільність); при дослідженні кори Землі за допомогою даних з кернів свердловинного 
буріння виникає задача відновлення внутрішньої структури між свердловинами. При 
цьому очевидним є той факт, що щільність грунту в різних точках кори є неоднорід- 
ною і найчастіше має розриви першого роду при переході від однієї складової кори 
до іншої (чорнозем, пісок, глина, граніт тощо). 

Весь розвиток обчислювальної та прикладної математики говорить про те, що 
використання кожної додаткової інформації про досліджуваний об'єкт може при- 
вести до більш точного 1 якісного відновлення цього об'єкта. Наприклад, в роботі |3| 
пропонується використовувати рівняння поверхні черепа людини і, таким чином, 
більш точно відновлювати внутрішню структуру тіла. Тому актуальною є розробка та 
дослідження математичних моделей розривних процесів. 


Постановка задачі 

Нехай досліджуваний процес описується функцією однієї змінної /(х) на інтер- 
валі |а,Б| та відомі точки розриву процесу ху, К -іп. Припускаємо, що хоча б в 
одному вузлі х; функція має розриви першого роду. Задані вузли розбивають ін- 


тервал (а, | на п частин. 


Метою даної роботи є побудова та дослідження математичної моделі заданого 
розривного процесу із заданими можливими точками розриву у вигляді розривного 
лінійного інтерполяційного сплайна. 


152 «Искусственньй интеллект» 272011 


Математичне моделювання процесів, які мають розриви... 2 


Побудова розривного інтерполяційного сплайна 


Визначення. Будемо називати розривним інтерполяційним лінійним сплайном 


на відрізку | ху, хо, ЮК -1,п-1 наступну функцію: 


ХХ 


509) 9 Зруб 4) А НН Аа ,Каїл-і, () 
ХХІ Хе 7 Ук 
де 
А; З Лу 0), Аа ЗЛ (ха 0). 
Теорема 1. Функція 5(х) - 5рі (х), К 1,п задовольняє наступним властивостям: 
Уру (ху - 0) з Аг, Уру- Оу - 0) - Ар: (2) 
Доведення. Слід врахувати, що зліва від вузла х; сплайн 5(х) задається формулою: 
5(2) щ Ура, 4) - А До ре о о 
Хр Х і з зо 
а справа - 
хх - хе 
500) є 5ру(о Аз я 
Хр Хр Храї Я 


Перевіримо виконання інтерполяційних умов (2). 


з  ХроХ - Хр ОХ Ж 
Уррл (ху - 0,4) 5 Ак НА -А, 
Хр Хр Хр 
з ХеОХро о до 0 Хе ХА з 
Урк (ху 1 0, А) я Арон ВА 
Хр Хі Хр ХЕ 


Теорема І доведена. 
Теорема 2. Якщо на кожному з інтервалів | ху, ху І, К - 1,п - 1 невідомі параметри 


А, Аа знаходити з умови 


тах тах |/(0)-5ру0)|-» тіп, (3) 
І«Ккеп- 1 хрохохр А 


то отримаємо розривний апроксимаційний сплайн, який є сплайном найкращого на- 
ближення. 

Доведення випливає з того, що кожний з елементів, який треба мінімізувати, 
дорівнює максимальному відхиленню наближуючого сплайна від функції /(х). Тому, 


при знаходженні параметрів з умови (2), отримаємо А, |, 4; з, які забезпечують най- 


менше відхилення. 
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Теорема 3. Якщо наближувана функція /(х) є розривною кусково-лінійною функ- 
цією з точками розриву хх, К -1,п і наближуємо її кусково-лінійним розривним 
сплайном 5р(х, 4), що визначається формулами (1), і невідомі Аг, А, | знаходимо з 


умови (3), то отримаємо точно наближувану функцію, тобто 
ур, 4) - /(9. 


Доведення. Нехай функція /(х) на Кк -му інтервалі має вигляд 


ХХ 


ловуава - НОВ» з каеінп-і. 


хуооХо хро 


Розглянемо максимум різниці функції /(х) на наближуваного сплайна (1) на К -му 


інтервалі: 
тах | Г(2)- 5ру С», 4) - 
ххьхні 
- о тах ФГОх -- 0) - 5Уру (ху -0, 4), (ху - 0) - Зрр(хру -0,40)) аз 
хе хр 


9 


ія 


- тах Ви о Ах 
хб6хроені 


В бу 


Знайдемо максимум за всіма інтервалами та мінімум отриманого максимуму 


Вб й : 


2 


пипо тах тах В У 
І«КкапоЇхрохохуцр | й 


Звідси випливає, що 


Ви у Ах -0, 


Вк Ак 2-05 874 84. 


Теорема доведена. 
Визначимо вигляд похибки наближення розривним сплайном (1) та оцінку 
наближення розривної функції побудованим сплайном, які наведені в роботі |4). 


Теорема 4. Якщо /(х) є С"|а,Б), г «1,2, то залишок ВК/(х) є /(2)- 509) на 
кожному інтервалі розбиття буде мати вигляд: 


моз | /ПДСоає, хеЇхьхьн 


ХЕ 


хохіи (07 
Ра -Ту 
сег хро (7-1) 
бохохр (Ху - 5) 
хуагхоо (7-1)! 


Р ХрОбЯХ ХХ, 


паза. 
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Теорема 5. Оцінка похибки наближення двічі диференційованої функції /(х) 
побудованим розривним інтерполяційним сплайном 5(х) - Уру (х) на кожному інтер- 


валі | хрухуці, К -1,п-1 має вигляд: 


ЛодеСхрхунів озир |/09-5руС9|« ча зно |/ 62 Ро 


хрехеху 


- 2 
Год є Срок) ЇКООЇ, є ЕН круто у 


Хені 


Зауваження. Якщо А, «А, з /(х), КеІ,п-і, то побудований розривний 


інтерполяційний сплайн вигляду (1) є неперервним лінійним інтерполяційним 
сплайном. 


Побудований розривний інтерполяційний сплайн 5(Х) у вигляді формули (1) є 
математичною моделлю розривного процесу, який описується функцією /(Х) з 
розривами у вузлах ху, К -1Іп. 


Приклад 1. Нехай задана функція /(х) на інтервалі |-1,|| з однією точкою 
розриву х 0) першого роду (рис. 1): 


задо г ОВ 


І, О0«хаї 


їСх) 0 


п ДВ 
| 0 1 
-1 х 1 
Рисунок І - Графічний вигляд наближуваної функції в прикладі 1 
Обираємо вузли сплайна: х| - - І, х, - 0, хз 


ронні значення функції у вузлах: 


-1. Вважаємо заданими односто- 


А «(а т 0)-0, 45 а Го є 0) 2-1, 
4, З 05 - 0) -1, 4 З (3 - 0) -0. 
Наближуємо функцію /(х) інтерполяційним сплайном вигляду: 
о о р 
хто хо -Х хі  -Іієдеб 
5р(х) з з і (9) 
р ХУ -хох» дені О«хаїЇ 
А5 4 , охусХаХЗ 
оо ЗС1Я0 


«Штучний інтелект» 272011 135 


л Литвин О.М., Першина Ю.І. 


Як бачимо, інтерполяційний сплайн точно наближує задану лінійну функцію, що і 
підтверджує вищевикладену теорію. 

Тепер побудуємо апроксимаційний сплайн у вигляді формули (1). Коефіцієнти 
матриці 4 знаходимо, застосовуючи теорему 2, тобто розв'язуємо мінімізаційну 
задачу: 


тах | /(х)- 5ру (х, 4) -» тіп. 
І«Ккхп-і1 4 


В результаті отримаємо таку матрицю А: 


47 Я) 


Тобто найкраще наближення заданої функції в прикладі І буде мати вигляд: 


хі, -І1сх«0 


Ур(х) - | 


хх, дах! 


що повністю збігається з побудованим інтерлінаційним сплайном. 
Приклад 2. Нехай задана функція /(х) на інтервалі |-1,1| з двома точками 
розриву х - - 0.5, х - 0.5 першого роду (рис. 1) 


хорі  ззфчже-ї8, 
ЛО0)г-хх, -0.5 «хх «0.5,. 
х2-1  05«ххї, 
1,96 
Ях) 
0,748 
-1 0 І 
0,98 х 0,98 


Рисунок 2 - Графічний вигляд наближуваної функції в прикладі 2 


Обираємо вузли сплайна: х; - - І, х» з - 0.5, хз з 0.5, хд -1. Вважаємо заданими 
односторонні значення функції у вузлах: 


А а Га 0)-2, Ж а 103 - 0) 0,5, 
45 а Го - 0) -1.25, Ж а ГО 0) 2-05, 
А з ГО 0) 2-05, А а Гц 0) 0. 
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Наближуючий інтерполяційний сплайн, згідно з формулою (1), буде мати вигляд 


(рис. 3): 
-1, Ух -- 0,5; -І1«х«-0,5; 
5р(х) зх; -0,5«х«0,5). 
1,5х - 1,5; 0,5«х«ї. 


ї(х) їх) 
-0А 
5(3, х, М/3) 


15. 
5(,х, М) 


І І 
-Р -09 -08 -07 -0.6 
-ї х -0.5 


а) 
Рисунок 3 - Графічний вигляд функції /(х) та наближуючого інтерполяційного 
сплайна Ур(х) на інтервалі: а) (-1; -0,5), б) (-0,5;0,5), в) (0,51) 


Максимальне відхилення функції /(х) від побудованого інтерполяційного сплайна 
дорівнює 
тах ГО - уро) я» 0,06. 
Тепер побудуємо апроксимаційний сплайн у вигляді формули (4), де коефіцієнти 
матриці А знаходяться з умови (3), тобто сплайн набуває вигляду (рис. 5): 
-1,5х- 0,46;  -І«х«-д,5; 
Ур(х) іх; - 0,5 « х « 0,5; 
1,5х -1,54; 0 5«х«Ї. 


1.999 


І І І І 
-1 -09 -08 -07 -06 
й Бе -05 


а) 
Рисунок 4 - Графічний вигляд функції /(х) та наближуючого апроксимаційного 
сплайна Ур(х) на інтервалі: а) (-1; -0,5;), б) (-0,5;0, 5), в) (0,541) 
Максимальне відхилення функції /(х) від побудованого інтерполяційного сплайна 
дорівнює 
тах |/О9 - уро) ях 0,02. 
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Висновки 


Таким чином, в роботі запропонована математична модель розривного процесу, 
що описується функцією однієї змінної з можливими розривами першого роду в зада- 
них вузлах, за допомогою розривного лінійного інтерполяційного та апроксимаційного 
сплайнів. 

Побудовану математичну модель можна буде використати для відновлення 
внутрішньої структури об'єктів, що мають різну щільність, в медичних, геологічних, 
космічних та інших дослідженнях. 
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О.Н. Литвин, Ю.Й. Першина 

Математическое моделирование процессов, имеющих разрьгвьг, 

с помощью разрьвньх интерполяционньгх сплайнов 

В работе предложен метод построения разрьтвного интерполяционного линейного сплайна для приближе- 
ния функции одной переменной, имеющей хотя бьг один разрьгв первого рода в узлах разбиения области 
определения функции, используя метод минимакса. Определен общий вид погрешности приближения 
функции построенной разрьтвной конструкцией в интегральном виде, и приведеньт оценки погрешности 
приближения в каждом интервале разбиения. Предложенньй метод можно использовать для математи- 
ческого моделирования разрьїювньтх процессов в медицинских, геологических, космических и других 
исследованиях. 
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